Long Time Tail of the Velocity Autocorrelation Function in a Two-Dimensional 

Moderately Dense Hard Disk Fluid 
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Alder and Wainwright discovered the slow power decay ~ t~''^^ (ci:dimension) of the velocity 
autocorrelation function in moderately dense hard sphere fluids using the event-driven molecular 
dynamics simulations. In the two-dimensional case, the diffusion coefficient derived using the time 
correlation expression in linear response theory shows logarithmic divergence, which is called the "2D 
long-time-tail problem". We revisited this problem to perform a large-scale, long-time simulation 
with one million hard disks using a modern efficient algorithm and found that the decay of the long 
tail in moderately dense fluids is slightly faster than the power decay (~ 1/t). We also compared 
our numerical data with the prediction of the self-consistent mode-coupling theory in the long time 
limit (~ l/{tVlnt)). 
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More than thirty-five years ago, Alder and Wainwright 
discovered the slow power decay ^ t^'^/^ (dtdimension) 
of the velocity autocorrelation function (VACF) in mod- 
erately dense hard sphere fluids using an event-driven 
molecular dynamics (EDMD) simulation [ij. In the two- 
dimensional (2D) case {d — 2), the transport coefhcient 
(i.e., the diffusion constant D) derived using the time- 
correlation function in linear response theory, the so- 
called Green-Kubo expression [3], shows logarithmic di- 
vergence. This means that conventional hydrodynamics 
do not exist, and this is called the "2D long-time-tail 
problem" [1, This discovery has greatly influenced 
the development of non-equilibrium statistical physics 
and liquid states; several theories based on a kinetic ap- 
proach [5|] and mode coupling theory (MCT) 0, 0] have 
been constructed. From the numerical perspective, since 
a longitudinal sound wave propagates at the speed of 
sound Cs in a system with periodic boundary conditions 
(PBCs) , the maximum correlation time for obtaining the 
true VACF in a numerical simulation is limited by the 
system size, i.e., the particle number, N. Therefore, we 
must perform a large-scale simulation to explore the ex- 
act form of the long time tail, which is not yet fully un- 
derstood p. 

Alder and Wainwright [l| obtained the VACFs in a 
hard disk system with about one thousand hard disks 
and found that long power decay occurred instead of 
ordinary exponential decay. Erpenbeck and Wood also 
examined this system in a direct simulation with sev- 
eral thousand hard disks 9]. In the 1990s, Frenkel and 
Ernst, Hoef and Frenkel, Naitoh et al., and Lowe and 
Frenkel investigated a systematic large-scale simulation 
using a lattice gas cellular automaton model and con- 
ducted theoretical analyses [l^, [HI, • Since their 
system is discrete in space, time, and velocity, the long 
time tail obtained from very efficient simulations can be 



compared with the predictions of the self-consistent mode 
coupling theory (SCMCT) [13, [3 j in which they showed 
a clear evidence of a good agreement between them. It 
has been very difficult to investigate the long tail using 
a direct hard disk molecular dynamics simulation using 
an event-driven scheme for several reasons. First, since 
a sound wave propagates through PBCs and the true 
VACF is disturbed, a large system with many particles 
is needed. Second, to discuss the tail of a long time corre- 
lation, we must perform a long time simulation for each 
parameter. Finally, to investigate the functional form 
of the tails in the VACF accurately, we need to average 
numerous statistical samples (ensembles) of independent 
physical properties such as velocities. Recently, the de- 
velopment of the computer and sophisticated modern al- 
gorithms [H, [13, [3 have enabled us to perform massive 
molecular dynamic simulations of a hard sphere system. 

In this study, we revisit the "2D long-time-tail prob- 
lem" and perform a large-scale, long-time, statistically 
accurate, systematic EDMD simulation with one million 
hard disks using a fast modern algorithm ^1^] . We found 
that the decay of the VACF in moderately dense flu- 
ids is slightly faster than (~ l/t), which seems to agree 
with the prediction of the SCMCT in the long time limit 

(~ i/itVhTt)). 

The Green-Kubo expression for the diffusion coefficient 
D is described by 



D 



{vx{t)vx{0))dt 



(1) 
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where (t) is the velocity of the tagged particle at time 
t and (•) indicates ensemble averages. In conventional ki- 
netic theory, the time correlation function {vx{t)vx(0)) 
decays exponentially as {vx{t)vxiO)) exp(— t/re), 
where te is the relaxation time. However, Alder 
and Wainwright [J] discovered the long time tail as 

{vx{t)vx{0)) ~ {t/TE)~'^^^ ■ In the case d = 2, the 
long time tail derives the logarithmic divergence as D 
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In (oo) + const., in which the hydrodynamics would break 
down. 

Our system consists of about one miUion elastic 2D 
hard disks {N = 1024^) placed in a L^; x Ly square box 
with PBCs. The basic units in this system are mass to, 
disk diameter oq, and energy k^T. The time unit can 
be described as y/rrm^^JhsT . Initially, the simulation 
systems for each packing fraction v — NTT{aQ/2)'^ / L^Ly, 
are prepared as the equilibrium state in a sufficiently long 
preliminary run, in which the density is uniform and the 
disk velocities have a Maxwell-Boltzmann distribution. 
The system evolves through collisions, with up to 10^ col- 
lisions per particle, using a modern fast algorithm based 
on event-driven molecular dynamics . To obtain accu- 
rate VACFs, we average samples using the particle num- 
ber N and velocities in two directions {vx^Vy). In this 
system, there are several time scales; the mean free time 
to, and the typical relaxation time te at which the de- 
cay gradually changes from the exponential to the power 
form. The fast longitudinal sound wave that results from 
a particle collision propagates with a sound velocity of Cs 
through the "artificial" PBCs. This artificial effect ap- 
pears in the VACF around time tmax ^ {L/'i)/cs- If we 
want to investigate pure VACFs at the thermodynamic 
limit, it is important to evaluate tmax (i-e., the sound 
velocity Cs). In all of the VACF figures in this letter, 
the time is scaled using the mean free time to ^nd the 
VACFs are normalized using the square of the initial ve- 
locity {v(S)f). 

Here, we summarize the theoretical estimation of the 
mean free time to and the sound velocitycg in a hard 
disk fluid. We refer to 2D Enskog theory ]l9| In Eq.(89) 
of Ref . [igj , the mean free path I as a function of number 
density n is described by 



V = 0.45 



2V2naoY{y) 8V2yY{y) ' 



(2) 



where y = (7r/4)riag and Y{y) = (1 - {7/16)y)/{l - yf 
(Enskog factor). The non-dimensional mean free path 
^*(y)(= l{y)/o.o) becomes 



8V2yY{y) ' 



(3) 



Therefore, the non-dimensional mean free time tp can 
be estimated easily using l*{y) and the mean particle 
velocity. The non-dimensional sound velocity c* in units 
of (oq, to, UbT) in a hard disk fluid for each density is also 
estimated using Enskog theory, and can be described as 



if + P + y 



dy' 



(4) 



where f{y) = l + 2yY{y) and df/dy = 2{\+y 18) / {\-yf . 
We confirmed that these theoretical values are in quite 
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FIG. 1: The VACFs for a typical packing fraction in terms 
of inverse time are shown. The particle number and packing 
fractions are fixed at iV = 1024 and v = 0.18,0.30,0.45, re- 
spectively. With these parameters, this figure corresponds to 
Fig. 3 in Ref. [ij. 



good agreement with numerical simulations for a wide 
range of densities, except around the solid-fluid transition 
density (the so-called Alder transition v ~ 0.70). 

Figure 1 shows the VACFs in a system with 32^ (= 
1024) hard disks, which is almost the same size as Fig. 3 
of Ref. [1], in which the maximum system size was 
N = 986. The packing fractions v (= 0.45,0.30,0.18) 
in Fig. 1 correspond to A/Aq {— 2, 3, 5) in Ref. fT], where 
A/Aq means the area A of the system divided by the 
area at the close packing density. (The relation be- 
tween v and A/Aq is V = 7r/(2V3(yl/Ao)).) We found 
that the curves obtained using our numerical simulation 
are quite smooth and consistent with Fig. 3 in Ref. jl|. 
In this small system, it is difficult to discuss the exact 
functional form of the tail and its exponents, since the 
region between te and tmax is quite narrow. 

To clarify the effect of sound wave propagation due to 
PBCs, we investigated the system-size dependence of the 
VACF by changing the particle number systematically. 
Figure 2 shows VACFs of various sizes at a fixed packing 
fraction [v ~ 0.30). We found clear evidence of artifi- 
cial disturbance of the sound wave by PBCs in VACFs, 
in which they deviate from the pure VACF around time 
tmax- For instance, in the case N = 16^ (= 256), it de- 
viates from the pure VACF at around t/to ^ 10, which 
corresponds to tmax for the system N = 16^ (= 256). 
The same results were obtained at other packing frac- 
tions. These facts confirm the validity of the sound ve- 
locity Cs estimated using Enskog theory in a hard disk 
fluid. The maximum size of our system is about one mil- 
lion particles {N — 1024^), for which tmax is estimated 
as t/to 698 a,t V — 0.30. Our results can be used to 
discuss the pure VACF of the thermodynamic limit at 
more than t/to 500. In the following, we analyse only 
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FIG. 2: The system-size dependence of the VACFs on the 
packing fraction f — 0.30. 



N=102r 



< 
> 




FIG. 3: The packing fraction dependence of VACFs from a 
dilute to a moderately dense fluid with a one-million hard disk 
system is shown. The vertical axis of the VACFs is multiplied 
by t/to in the semi- log plot. 



the data of the VACFs with the maximum system size 
{N ^ 10242 = 1048576). 

Under the assumption that the decay of the long tail 
has the power form, as indicated by Ref. [l|, we first 
fitted the numerical data between te and tmax with a 
function in the form ~ (t/to)" by changing the parameter 
a. We found that the exponent a obviously deviated 
from —1, and had values of a ~ —1.08, —1.06, and —1.03 
for I' = 0.18, 0.30, and 0.45 respectively. The numerical 
data of the long time tails in a moderately dense system 
obtained from our simulation seem to decay faster than 
the power form 1/^)- 

To investigate the functional form of the long time tail 
in detail, we show the semi-log plot of the VACFs for the 
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FIG. 4: The packing fraction dependence of the VACFs from 
a moderately dense fluid to a solid with a one-million hard 
disk system is shown. The vertical axis of the VACFs is mul- 
tiplied by f/to in the semi- log plot. 



packing fraction in Figs. 3 and 4. The vertical axes are 
multiplied by tjto to show the deviation from the con- 
ventional prediction of decay (tjto)"'^ more clearly. In 
the stage of exponential decay < te)-, the universal 
curves of the VACFs for each packing fraction have a 
maximum around tjto ^ \. In the case of a dilute gas 
(y = 0.05), the tail of the VACF (i/io > 10) seems close 
to a flat line (the power decay with a ^ —1). By con- 
trast, in a moderately dense gas [y = 0.15 ~ 0.35), the 
decay of the VACFs is faster than the conventional power 
decay. In Fig. 4, the "back scattering" effect gradually 
becomes more dominant than v 0.45 around time te- 
At = 0.45 ^ 0.60, this tendency becomes more remark- 
able. In a dense fluid at ^ = 0.65, VACF takes a negative 
value. At ^ = 0.70 ~ 0.75 (more than the Alder transi- 
tion point), the effect of the solid- fluid transition causes 
a drastic change in the VACF, even when the decay is 
exponential. We found that these results reveal new nu- 
merical results for the long time tail that show a peculiar 
density dependence. These facts have never been inves- 
tigated in previous works [l|, since the true VACF can 
only be seen within the time scale t/t^ 30. 

From a theoretical perspective, another possibility for 
decay faster than the power form (~ (1/^)) has al- 
ready been proposed. This is weak logarithmic diver- 
gence based on the self-consistent mode coupliirg theory 
(SCMCT) [li,[ll. Ernst et al. @ predicted the VACF 
Coit) based on MCT in a 2D hard disk fluid with the 
shear viscosity as 



Coit) 



8m-K{v^is + D) 



(5) 



However, this solution results in a contradictory conclu- 
sion in assuming a finite value of the diffusion constant 
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FIG. 5: The VACFs in a moderately dense hard disk fluid are 
compared between numerical simulations and the theoretical 
predictions using the Simple MCT @] (broken line) and self- 
consistent MCT [13,111 (solid line). 



D as D diverges in the framework of linear response re- 
lations (Green-Kubo expression f3|). To avoid this prob- 
lem, Kawasaki [3l and Wainwright et al. [15| proposed 
the SCMCT independently, which is described by 



Coit) 



ksT 
16mn 



(6) 



In this case, the diffusion coefficient D{t) diverges weakly 
at the long time limit as D (In (oo) -I- const.) ^ . Wain- 
wright et al. [l5|] stated that the two decays in Eqs. (5) 
and (6) cannot be distinguished by their simulation 
within the maximum correlation time t/to ^ 30. Figure 
5 compares numerical simulations of VACFs in a mod- 
erately dense fluid and the theoretical prediction of the 

SCMCT as a function in the form (t^/\n (t)^ ) with 



a fitting parameter /3. Although the decay can be dis- 
cussed within t/tf) ^ 10^, the two functional forms at the 
long time limit seems to converge. 

To summarize, we revisited the "2D long-time-tail 
problem" using a modern fast event-driven molecular dy- 
namics simulation [l^. We completely reproduced the 
results of Alder and Wainwright, who discovered the 
power form of the long time tail [li]. Compared with 
previous works, we essentially differ in having explored 
in detail a larger sample, longer correlation time, and 
more accurate VACFs. We found that the tail of VACFs 
seem to decay as the power form (~ t~^) at the low 
density {v — 0.05). The remarkable results obtained 
from our simulation are the discovery of the decay of 
VACFs in a moderately dense hard disk fluid, which is 
faster than the conventional power decay. We also com- 
pare the prediction of SCMCT [H, [H seems consistent 
with our numerical results at a long time limit exceeding 
t/to ~ 100. Although the coefficients of the prediction in 
the SCMCT change with the fitting parameter in Fig. 5, 
the tails seems to converge in the long time limit. We 
conclude that a simple description of the decay of the 
power form over a long time limit might not be correct, 
at least in a moderately dense hard disk fluid. Based 
on our results, further studies should reconsider 2D hard 
disk fluid using both extensive numerical simulations and 
the theoretical derivation of the transport coefficient (dif- 
fusion, viscosity, and heat conductivity) based on the lin- 
ear response, kinetic theory, mode-coupling theory, and 
hydrodynamics. 
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